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Abstract. We present preliminary results of our work on blazar unifi- 
cation. We assume that all blazars have a broad line region (BLR) and 
that the size of the BLR scales with the power of the source in a manner 
similar to that derived through reverberation mapping in radio quiet ac- 
tive galactic nuclei (AGNs). Using a self-consistent emission model that 
includes particle acceleration we show that according to this scaling, in 
weak sources like MKN 421, the inverse Compton (IC) scattering losses 
are dominated by synchrotron-self Compton scattering (SSC), while in 
powerful sources, like 3C 279, they are dominated by external Compton 
(EC) scattering of BLR photons. In agreement with other workers, we 
show that even in the powerful sources that are dominated by EC scat- 
tering, the hard X-ray emission is due to SSC. Finally, we show that this 
scaling reproduces well the observed sequence of blazar properties with 
luminosity. 



1. Introduction 

Blazars have been shown to exhibit a sequence of properties as a function of 
source power (Fossati et al. 1998). As the source power increases, the emission 
line luminosity and the ratio of Compton to synchrotron luminosity increase, 
while the synchrotron peak frequency v s and the IC peak frequency decrease. 
Recent multiwavelength studies support this scheme (e.g. Kubo et al. 1998), 
although the result that the ratio of the Compton to synchrotron luminosity 
increases with source luminosity suffers from limited statistics, and should be 
only considered tentative. 

The initial division of blazars into flat spectrum radio quasars (FSRQs) and 
BL Lacertae objects (BLs) was based on the equivalent width (EW) of the broad 
emission lines. Sources with EW > 5 A were classified as FSRQs and sources 
with EW < 5 A as BLs. This difference in the EW of the emission lines has 
been interpreted as absence of a substantial BLR in BLs, and has been used 
to advance the idea that in BLs the GeV-TeV emission is due to SSC, while 
in FSRQs the GeV emission is due to EC. However, the EW criterion does not 



1 



2 



Georganopoulos, Kirk, and Mastichiadis 



correspond to a dichotomy (Scarpa & Falomo 1997), because a weak BLR is 
present in BLs as well. Additionally, it is not the BLR luminosity Lb that is 
relevant for the EC luminosity, but the BLR photon energy density Ue that is 
measured in the comoving frame of the non-thermal source. One needs to know 
not only Lb, but also the BLR radius Rb and the location of the non-thermal 
emitter relative to the BLR in order to quantify the relative importance of SSC 
versus EC emission. 

The size of the BLR has been measured only for radio quiet AGNs through 
reverberation mapping (Kaspi et al. 2000), and it has been found to scale 
with the ionizing luminosity L acc of the acrretion disk as Rb oc LjJcc- ^ was 
additionally found that, assuming Keplerian motions for the BLR clouds, the 
mass of the central object scales with the ionizing luminosity as M oc l" c 5 c . 

Here we assume the scalings apply also to radio loud AGN and examine the 
implications for blazars. Following the formalism of internal shocks propagating 
in a conical jet (Rees 1978), we scale the location Duob and the radius Ruob of 
the non-thermal spherical emitter with the mass M of the central object. We 
assume that Lb is a fraction of L acc and scale the kinetic luminosity and the 
Poynting flux of the blob with L acc . 



2. The scaling and its application 

Under this scheme, the radius of the BLR scales as Rb oc L^J c , and the distance 
of the blob from the center of the system scales as Duob oc M oc L®f c . Therefore, 
we assume that for powerful sources the blob radiates from inside the BLR, 
Dbiob < Rb- However, if we consider increasingly weaker sources, Rb falls 
faster than Duob and, below some critical luminosity, the blob radiates from 
outside the BLR. For a blob inside the BLR, Ue oc (Lb/R b )^ 2 , where T is 
the bulk motion Lorents factor of the blob. For a blob outside the BLR, the 
solid angle subtended by the BLR, as seen by the blob, is reduced. This affects 
Ue in two ways. The first is the usual geometric 1/r 2 factor, and the second 
is a relativistic de-boosting that, for 3> Rb-, gives Ue oc (Lb/ D^ lob )/T 2 

(Dermer & Schlickeiser 1994). The second effect dominates if the blob is located 
less than a few Rb away from the BLR. For example, for T = 10, Ue at a 
distance of 10 Rb is smaller than Ue inside the BLR by a factor of 10 2 due to 
the usual 1/r 2 and by a factor of ~ 10 4 due to relativistic de-boosting. 

We apply this scaling using a model (Georganopoulos & Kirk 2000) for 
the blob which includes both particle acceleration and radiative losses with IC 
scattering losses treated in the Thomson regime. Particles are accelerated in an 
acceleration zone and escape into a radiation zone, before they eventually escape 
out of the system. 

This simplified picture is intended to represent particles which are acceler- 
ated by a shock front, escape it, and radiate downstream before the compressed 
plasma re-expands. Whilst undergoing acceleration, particles simultaneously 
suffer synchrotron losses and losses by Compton scattering of both photons of 
external origin (EC) and synchrotron photons (SSC) from the radiation zone. 
In the radiation zone, EC, SSC and synchrotron losses occur until the particle 
finally escapes after a time t esc = 3Ru b/c- The electron energy distribution 
(EED) is computed self-consistently, taking account of these processes. 
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Figure 1. Sizes (lower panel), energy densities (middle panel), and 
the break of the EED as a function of the accretion disk luminosity. 



Our assumption about the scaling of the Poynting flux implies a constant 
magnetic field, which we take to be B = 0.2 G. The blob is assumed to move 
with a Lorentz factor T = 15 at an angle 9 = 1/T with respect to the line of 
sight. Low energy particles (70 = 10) are injected into the blob at a rate of 
Q = 0.1 cm -3 s _1 . Also Lb = 0.01L acc . For L acc = 10 44 erg s _1 , reverberation 
mapping gives Rb ~ 8.6 10 16 cm. We assume that at this luminosity Rb = L>bi b 
and R b iob = 2 10 15 cm. 

In Fig.l, we let L acc vary from 10 erg s _1 to 10 46 erg s _1 , and follow 
the behavior of the system. In the lower panel Rb, Dbiobi an d Rbiob are plotted 
as a function of L acc . For weak sources, the blob is located outside the BLR 
{Dbiob > Rb)- As the source power increases, the blob gradually approaches and 
eventually enters the BLR. In the middle panel of Fig.l we see how this affects 
Ue (the magnetic field energy density Ub is constant under the adopted scaling). 
For weak sources, Ue is much smaller than both Ub and the self-consistently 
derived synchrotron photon density Us- Gradually Ue increases, and eventually, 
for bright sources, dominates over both Ub and Us- In the upper panel of Fig.l 
we plot the self-consistently derived break jbreak m the EED. Note that ^break 
decreases as the power of the source decreases. 

In Fig. 2 we plot basic observable quantities as a function of L acc . In the 
lower panel we plot the synchrotron luminosity L5, the SSC luminosity Lssc-, 
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Figure 2. Luminosities (lower panel), peak frequencies (middle 
panel) and the ratio or synchrotron to accretion disk luminosity and 
IC to synchrotron luminosity (upper panel). 



and the EC luminosity Eec- At low powers, Lec is much weaker than L$sc 
and Ls that are roughly equal. As the source power increases, Lec gradually 
dominates over L$sc an d Ls, and we end up with an EC dominated source. In 
the middle panel we plot the peak frequencies of the three emission components. 
Note how the synchrotron peak frequency us decreases as the source power 
increases. Similar behavior is also seen for the SSC peak frequency fssc and 
the EC peak frequency vec- Finally in the upper panel of Fig.l we plot the 
ratio Ls/L acc (solid line) of the synchrotron to accretion disk luminosity and 
the ratio Lc / Ls of the IC luminosity to synchrotron luminosity. Ls dominates 
over L acc in weak sources in agreement with the lack of a thermal component 
and weak emission lines in weak sources like BLs. L acc becomes more significant 
for powerful sources, again in agreement with the strong emission lines of FSRQs 
and the accretion disk signature observed in some FSRQs (e.g. in 3C 279, Pian 
et al 1999). The dominance of the Compton over the synchrotron emission 
increases as the source power increases, in agreement with observations (e.g. 
Kubo et al. 1998). Note, however, that the observational case for an increasing 
Compton dominance with source power is based on poor statistics due to the 
limited sensitivity of the CGRO. 
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Figure 3. Synchrotron luminosity versus synchrotron peak frequency 
for a set of blazars (Sambrunaet al. 1996, Kuboet al. 1998). Overlayed 
are the the model tracks as a function of L acc for a range of observing 
angles. 

In Fig.3 we plot the synchrotron luminosity L s versus the synchrotron peak 
frequency u s for the blazars studied by Sambruna, Maraschi & Urry (1996) and 
Kubo et al. (1998). On top of the data points we plot the model tracks as a 
function of L acc for a range of observing angles. Note that these tracks represent 
two model derived quantities, Ls and vs as L acc varies. The model covers rather 
well the observed parameter space. Given the track of the model under an angle 
9 = 0°, we do not expect that any powerful (Ls ~ 10 47 erg s -1 ) sources with 
high peak frequencies vs ~ 10 17 Hz exist. The discovery of such sources would 
pose a serious problem for this model. 

3. The luminosity scaling 

In Fig. 4 we plot the peak luminosities of the three emission mechanisms for 5 
model sources with L acc ranging from 10 42 up to 10 46 erg s -1 . The behavior 
of the model is in good agreement with the general characteristics of the ob- 
served luminosity sequence (compare with Fig. 12 of Fossati et al. 1998). The 
synchrotron peak frequency decreases from 10 18 Hz down to 10 12 Hz, as the 
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Figure 4. The model luminosity sequence. The synchrotron (circles), 
SSC (triangles), and EC (squares) luminosity as a function of frequency 
for 5 model sources of different intrinsic power. The broken lines link 
points from a single model. 



synchrotron power increases. The Compton peak frequency in weak sources is 
in the TeV regime (~ 10 25 ~ 26 Hz) and it is due to SSC, while in powerful sources 
it is in the GeV regime (« 10 23 Hz) and it is due to EC. According to the model, 
although in these bright sources the energy losses are dominated by EC, the hard 
X-ray emission is due to SSC, in agreement with resent observations (Kubo et 
al. 1998). 

We note here that Fig.4 corresponds to sources that are oriented at an angle 
9 = 1/r. In reality one should expect a range of angles, which will give rise to 
significant scattering around the presented trend. This scattering is visible in 
the work of Fossati et al. (1998) and indicates that a proper unification scheme 
for blazars should include the effects of orientation. 
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